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ABSTRACT 

We study a large set of high spatial resolution optical rotation curves of galaxies with 
the goal of determining the model parameters for a disk embedded within a cold dark 
matter (CDM) halo that we model either with a Navarro, Frenk & White (NFW) 
profile or pseudo-isothermal profile. We show that parameter degeneracies present in 
lower resolution data are lifted at these higher resolutions. 34% of the galaxies do not 
have a meaningful fit when using the NFW profile and 32% when using the pseu- 
doisothermal profile, however only 14% do not have a meaningful fit in either model. 
In both models we find correlations between the disk baryon fraction fd and the spin 
parameter of the halo A', between fd and the dark halo mass M 2 oo, and between M200 
and the concentration parameter c. We show that the distribution of the concentration 
parameters c, for a NFW halo, is in good agreement with CDM predictions; no sig- 
nificant galaxy population is found with very low values of c. The overall distribution 
of A' is in good agreement with theoretical predictions from hierarchical tidal torque 
theory. The whole sample is also well fitted by a pseudo-isothermal dark halo with a 
core, but the size of the core is rather small (6% of the virial radius or smaller; for 
70% of the sample the core size is less than 2 kpc) . Thus we conclude that the profile 
of dark matter is steep (r" 1 or steeper) down to this radius; large dark matter cores 
(and therefore very low dark matter central densities) seem to be excluded. LSBs tend 
to have higher values of A' for a given fd and lower values of c for a given mass than 
HSBs. In an appendix we give some useful formula for pseudo-isothermal profile halos 
and discuss in detail the issue of parameter degeneracies. 

Key words: cosmology: theory — galaxies: formation — galaxies: spiral — galaxies: 
kinematics and dynamics 



1 INTRODUCTION 



density cusp predicted by CDM (e.g., McGaugh fe dc Blok 



(1998j ); |Moore et al. (1999| ); [Dalcanton fc Bernstein (200(1 )) 



The Cold Dark Matter (CDM) paradigm for structure for- 
mation has proved remarkably successful in explaining the 
observed large-scale properties nf the universe, suc h as the 



Furthermore, numerical SPH simulations of disk forma- 
tion within dark halos fail to match the zero-point of the 
observed Tully-Fisher. One possible solution is t o lower the 



central concentra tion of the dark matter halos ( Navarro 



abundajqge and clustering of rala.xics and clusters ( Peacock Steinmetz 2000b[) 



et al 2001: Verde et al. 2001 



Lahav ct al. 2001 ), the sta- 



Phillips et al 



These problems have stimulated numerous proposed as- 
trophysical solutions, as well as mo difications of the fun- 



tistica l properties of the Lya forest (e.g 

(2001 )) and the power spectr um of the cosmi c microwave 

ha,cWr{nnd anisotronies (e.r.. [laffe et al (Tmfa ). However. (2000l): iKamionkowski fc Liddle (20001): IColombi. Dodclson 



damen t al CDM paradigm itself (e.g., I Spergel fc Stcinliardt 



a numt or of nuzzling- discrepancies remain when CDM nre- 



fc Widrow (19961): Bode, Ostrikcr fc Turok (200l|): iGoodman 



(20001 ); |Hu, Barkana fc Gruzinov (2000[ ); |Ccn (200 1| )) 



diet ion j are extrapolated to small scales. Chief among them 
are: (i) The substructure problem. CDM over-predicts the 
number of satellites around a Milky- Way sized galaxy by a n 



order of magnitude ( Klypin et al. 199£ ; Moore et al. 1999 ) . 
(ii) The density profile problem. Observed rotation-curves of 
dwarf and lower surface brightness galaxies suggest that the 
inner regions have a constant density core rather than the 



In this paper we once again examine the density profile 
problem by fitting CDM models to observed rotation curves 
of spiral galaxies. Our study differs from previous work in 
two respects. Firstly, we use optical rotation curves rather 
than the HI rotation curves us ed in previous studies which 
are plagued by beam smearing ( Swaters et al. 2000 ; van den 
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Bosch et al. 2000). Even if the effects of beam smearing 
are neglected, the relatively large errors and limited spatial 
sampling of HI rotation curves imply that they cannot be 
used t o discriminate between constant density cores and r _1 
cusps (van den Bosch & S waters 2001). By contrast, the op- 



tical rotation curves we use are free from beam smearing, 



halo and within a pseudo-isothermal profile halo (i.e., an 
isothermal sphere with a constant density core). The choice 
of a pseudo-isothermal halo is motivated by recent claims 
in the literature that the dark matter profile of low-s urface 
brightness galaxies (LSB, e.g., de Blok et al. (2001 )) and 
some h igh-surface brightness g al axies (HSB, e.g. , tBalucci 



have smaller errors and higher spatial resolution. We show (2001); Borriello & Salucci (2001); Salucci fc Burkert (2000[ )) 



that, with this superior data, various parameter degenera- 
cies present when fitting HI rotation curves can be lifted. 
This allows us to distinguish between core and cusp-like in- 
ner profiles. Secondly, while many studies have focused on 
relatively small samples of dwarfs (but see Navarro (199S ) ) , 
we use a large sample (400 galaxies) spanning a wide range 
in luminosity and surface brightness. 

Our goal is to determine the best fitting model pa- 
rameters for a disk within the CDM profile proposed by 
Navarro Frenk and White (1997; NFW) and within a 
pseudo-isothermal profile halo and study the distributions 
of the recovered disk parameters. 

We find that the NFW profile provides a good fit to 66% 
of the galaxies in the sample, with a distribution of recov- 
ered concentration parameters and spin parameters broadly 
consistent with that predicted by CDM numerical simula- 
tions. When the sample is fitted with an isothermal profile 
with a core, 68% of the galaxies are well fitted whithin this 
model and the best fit parameters favor cores with small 
sizes (below 6% of the dark halo virial radius). Many rota- 
tion curves that have no meaningful fit in one dark matter 
profile are well fitted by the other. There is no meaningful 
fit in either model for only for 14% of the sample. This is 
consistent with the inner dark matter profiles being steep 
(slope of -1 or steeper) down to radii that are few percent 
of the virial radius. 

In both models the recovered baryonic mass to light 
ratios are broadly in agreement with predictions from syn- 
thetic stellar populations. 

This paper is organized as follows. In section 2 we 
present the set of rotation curves we analyze and the two 
models we use to fit the data (exponential disk embedded 
in a dark matter halo with a NFW profile or a pseudo- 
isothermal profile). These two models have four free param- 
eters, but since for most of the galaxies the scale length of 
the disk is given, we report the analysis for four free pa- 
rameters in appendix II and present in the main text the 
results for the fit with three free parameters. In section 3 we 
illustrate and discuss the correlations we find between the 
best fit parameters for the two models. We discuss the gen- 
eral implications of these findings in section 4 where we also 
compare the derived baryonic mass-to-light ratios with the 
range allowed by stellar populations. Finally, we conclude 
and summarize our results in section 5. In appendix I we 
develop the expression for the rotation curve of an exponen- 
tial disk embedded in a pseudoisothermal halo. In appendix 
III we study the degeneracies among the parameters of the 
models. 



2 ROTATION CURVES 

We study a large set of observed rotation curves of galaxies 
with the aim of det ermining the model parameter s for a 
disk within a NFW ( |Navarro, Frenk fc White 1997j) profile 



is better represented by a flat core rather than the steeper 
profile found in CDM N-body simulations. 



2.1 The observational sample 

We use three sets of observed rotation curves. The first is 



the catalog compiled by Courteau (1997) (Courteau sample 



which consists of optical H a long-slit rotation curves for 300 
Sb-Sc UGC galaxies and i?-band photometry. The second is 
a sample of 74 s piral galaxies (Sa-Sd) observ ed again in the 



optical {H a ) by Palunas fc Williams (2000) (PW sample) 



For this sample I-band photometry and the velocity field 
in two dimensions is available. Finally, we also consider the 
recently published optical (H a ) rota tion curves for 26 Low 



Surface Bri ghtness (LSB) galaxies by McGaugh, Rubin fc de 



Blok (2001). These consist of long slit H a observations for 
26 galaxies and B-band photometry for a smaller sample. In 
total we have 400 optical rotation curves. 

All three samples provide inclination corrected curves 
with their corresponding errors. The Courteau and PW sam- 
ples have surface photometry and give values for the stellar 
disk scale lengths, but we do not have disk scale lengths for 
the LSB sample. 

The sample of galaxies studied here is somewhat het- 
erogeneous: we use LSB and HSB galaxies; some curves have 
much smaller error-bars than others, and the sampling range 
of the rotation curve varies significantly: in particular, some 
curves do not extend to large enough radii to show the flat- 
tening of the rotation curve. The diversity of the sample is 
both a strength as well as a weakness: it allows us to probe 
general statistical trends across a wide range of galaxy lu- 
minosities and surface brightness. 



2.2 The models 

Many authors have investigated gala c tosynthesis models 
(jDalcanton^^SjDcrgcl^j^Su 



Hawkins & Padoan 1997|; |MoJVlao & White 1998t Jimenez 



Padoan^ Matteucci & Heavens 1998]; Somerville & Primack 



199S|; |van den Bosch 2000|; |Avila- Reese fc Firmani 2000|; [Fir 



tnani & Avila- Reese 2000; Navarro & Steinmetz 2000b) in 



which the properties of disk galaxies are determined pri- 
marily by the mass, size, and angular momenta of the ha- 
los in which they form, and which may contain the effects 
of supernova feedback, adiabatic disk contraction, cooling, 
merging, and a variety of star-formation recipes. In this pa- 
per, we are mainly concerned with the dynamical properties 
of the disk/ halo system which affec t the rotation curve, i n 
the spirit of Dalcanton et al. (1997); Jimenez et al. (1997), 
and Mo, Mao fc White (1998[ ) (hereafter MMW). We model 



spiral galaxies as exponential disks (i.e., disks with surface 
density profiles S(r) = Eo exp(— r/Rd), where Rd is the 
scale length of the disk) embedded within CDM halos which 
we model either with NFW profiles or pseudo-isothermal 
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Figure 1. NFW model. The location of the sample of galaxies in different projections of th e three dimension al parameter space (Af200i 
fri. c) when fitted with a N FW dark halo model. The diamonds correspond to galaxies from CJourteau (1997) and the triangles are from 
Palunas & Williams (2000| ). The solid and dashed lines in the top-right panel are the N-body simulations predictions for the C-M200 
correlation (sec text fo r more details) Also plotted (bottom-right panel) is the deduced value for A' vs. fd- The solid line is the least square 
fit from 14 dwarfs by Burkert (200C] ), in good agreement with the correlation found here from 362 galaxies. Note also the correlation 
between fd and A/200- 



spheres. In both models we assume that the baryons initially 
follow the same density profile as the dark matter and then 
cool and settle into an exponential disk in a dynamical time. 
We assume that the baryons that do not settle in the disk 
follow the dark matter distribution. The disks have a mass 
fraction fd = Mdi a k/M2oo ^ f2f,/f2 m . Here Mdisk denotes the 
baryonic mass that settle into an exponential disk, M200 de- 
notes the mass within a radius -R200 that encloses an average 
density 200 times the critical density, il m denotes the den- 
sity parameter and fit the baryon contribution to the den- 
sity parameter. Throughout, we assume a ACDM cosmology 
given by: (fi m , fi A , Q b , h, a 8k -i ) = (0.3,0.7,0.039,0.7,1.0). 
In these models, the rotation velocity is given by: 



(1) 



where the bary onic disk rotation velocity is (Binney & 
Tremai jie (19871 ), p. 77) 



V![r) = ATvGZoRdy 2 \h{y)K (y) - h{y)K^y)} , 



(2) 



with y = r/(2R d ), E = f d M 200 /(2TTR 2 d ), I n and K n de- 
noting the modified Bessel functions of the first and second 
kind. The contribution to the rotation curve from the dark 
matter is: 



VcdmO) = GM CD M(r)/r, 



(3) 



where M(r) denotes the mass enclosed within the radius r. 
For a model to be fully specified, we require expressions for 
the disk scale length R d and for M(r). 

The NFW fit to the density profile of dark matter found 
in collisionless N-body numerical simulations is: 



p(r) 



: Pcrit 



(r/r e )(l+r/r c ) a 



(4) 



where p CT1 t is the critical density (p cr it = 277. 3/i 2 Mq 
kpc -3 ) and r c is roughly the break radius at which the 
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slope of the profile changes from — f to —3. The con- 
centration parameter c is defined as c = R.2oo/r c , where 
-R200 = (M2oo/(4/37r200p C rit)) 1/ ' 3 and M200 is assumed to 
be the total mass of the dark matter halo. The dark matter 
mass within an enclosed radius r is: 

— 1 +ln(l + cx) 



M(r) = 47r / o crit <W 

where x = r / R200 
80 and c is: 



1 



(•>) 



So = 



200 



From equation the relation between 



(6) 



3 ln(l + c) - c/(l + c) 

Disk formation affects the dark matter profile by the adia- 
batic contraction it induces in the inner regions of the dark 
matter halo. As in MMW, we model this by assuming that 
the halo responds adiabatically to the slow contraction of 
the disk, and therefore that the angular momentum of dark 
matter particles is conserved. This yields the relation be- 
tween the initial radius Ti and final radius r where a dark 
matter particle ends up: 



GM f (r)r = GM(r 



(7) 



where M(ri) is given by equation (g), and Mf(r) = 
E D exp(— r/Rd) + M(r;)(l — fd) is the total final mass within 
r (see MMW for details). 

Assuming angular momentum conservation, MMW de- 
veloped a fitting formula for Rd(c, M200, A', fd) which we use 
(see their equations (28) and (32)). Note that we use the 
modified spin parameter A' = (jd//d)A, where j d = Jd/J 
and Jd, J are the total angular momentum of the disk and 
halo respectively. If the specific angular momentum of the 
disk is equal to that of the halo, then jd = fd and A' = A. 

In the NFW model, for a given M200, the overdensity 
So and thus the concentration parameter c depends on the 
collapse redshift. There is therefore a correlation between c 
and M200 , since more massive halos collapse at lower redshift 
than less massive ones. However, there is a fai rly wide sca t- 
ter in the distribution of c for fixed M200 (e.g., Jing (200C)). 



Furthermore, there is growing observational evidence that 
the distribution of concentr ation parameters might not b e 



as predicted b y CDM (e.g., Navarro fc Steinmetz (2000b) 



Keeton (2001)); one of the goals of this paper is to explore 
this question further. We therefore keep c as a free parame- 
ter. Thus, there are four free parameters in our model: the 
spin of the dark halo A', the mass of the dark halo M200, 
the concentration parameter c and the fraction of the total 
mass in the disk fd. 

For the pseudo-isothermal sphere the profile of the dark 
matter is: 



P = Po 



1 + 



(0 



(8) 



where po denotes the finite central density. We have devel- 
oped a similar analytic model to compute the properties of 
an exponential disk in this profile (see Appendix I). In this 
case, we do not include the adiabatic compression of the 
halo due to the baryons, since we want to use a constant 
density inner core for the final dark matter distribution by 
fiat. In the pseudo- isothermal profile case, the CDM rotation 
velocity is: 



Vcdm = ^Gporl 



1 arctan 




(9) 



logioM 



Figure 2. NFW model. The distribution of c (left panel) and 
A' (right panel) deduced from NFW halos fits to rotation curves. 
Also over-plotted as continuous lines are the theoretical predic- 
tions for each of these parameters. Note that the derived values 
of c are broadly in agreement with the predictions from N-body 
simulations (see text for more details). 



A fitting formula for the disk scale length Rd as a function 
of A/200, c, fd and A' for this profile is given in Appendix I. 
Any rotation curve can be fitted with this profile using four 
parameters: A', M200, r c , f d . 

For each galaxy rotation curve we explore the whole 
likelihood surface and the best fitting model is then obtained 
using a standard \ minimization. 

We discuss in detail the values for the disk parameters 
found for these two models (NFW and pseudoisothermal) 
when using all four free parameters in Appendixes II and 
III. In this section we eliminate one of the four free param- 
eters (A') by using the observed scale lengths of the expo- 
nential disks Rd. Then the only remaining free parameter 
for the baryonic disk is fd. Since we do not have Rd for the 
LSB sample and for some of the galaxies in the PW sample, 
the analysis with three free parameters is performed on the 
Courteau and PW samples (in total 362 galaxies). The disk 
parameter distributions for the four parameters fit is quite 
similar to that for the three parameter fit. We discuss the 
locus of LSB galaxies in the parameter space in Appendix II 
and show that they follow the same trends as the Courteau 
and PW samples. 



3 BEST FIT MODELS 

3.1 Fits to disk models within NFW profiles 

Since we use the observed (stellar) scale lengths of the bary- 
onic disk (Rd), the free parameters in our model are: M200, 
fd and c. For each galaxy rotation curve we explore the en- 
tire likelihood surface within the following boundaries: 10 9 < 
M 2O o/M < 7 x 10 13 , O.K c < 100 and 0.001 < f d < 0.36. 
The choice of the boundaries is justified as follows: for M200 
by the typical velocities at large radii in the sample, for c 
by theoretical ansatzs given by N-body simulations and an- 
alytic predictions (e.g., MMW); and finally for f d a realistic 
upper bound can be given by nucleosynthesis constraints 
on the baryon fraction f d ~ 0.2 (e.g. laffe et al (2001 ); Efs- 
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Figure 3. Pseudo-isothermal model. Same as fig. |l| but when using the pseudo-isothermal dark matter halo. Note the similarity with 
fig- 



tathiou et al. (2001)); a realistic lower bound is more difficult 
to set, but it is reasonable to believe that at least few per- 
cent of the baryons in the halo cool to form the disk. We set 
a higher value fd = 0.36 and a lower value of /<j = 0.001 to 
explore parameter space more fully. The best fitting model is 
then obtained by using a standard inverse variance weighted 
X 2 minimization^. We compute the value of \ 2 on the entire 
grid of points to ensure that we converge on a global, rather 
than local, minimum in \ 2 ■ The discreteness of points in our 
solution reflects the finite spacing of our grid. 

Figure ^ shows the position of each galaxy for all the 
projections of the three dimensional parameter space (M200, 
fd and c). In addition, for the best fitting values of M200, fd 
and c we have computed A' and plot it against fd- 

For 122 galaxies out of the 362 the x 2 minimization pro- 
cedure does not find a global minimum within the bound- 
aries for the disk parameters. This means that for ~ 1/3 (i.e., 



The best fitting model to all 400 rotation curves can be found 



34%) of the galaxies in the sample the model with a NFW 
profile for the dark matter fails to reproduce the observed 
rot ation curve. This is in agreement with a previous study 
by van den Bosch & S waters (2001). They analyze (with 
a different procedure and different modelling assumptions 
from ours) 20 rotation curves of late type dwarfs and find 
that for 6 of these galaxies (i.e., 30%) no meaningful fit can 
be obtained. It is important to note that our analysis uses 
a different approach, a much larger sample of galaxies that 
is not limited to dwarfs, and still concludes that the same 
percentage of galaxies are well fitted by a CDM profile. 

There are some correlations among disk parameters 
worth commenting on. The top-left panel of fig. [l] shows 
an anti-correlation between M200 and fd, with a clear zone 
of avoidance for low masses and low spins (we discuss the 
signific ance and possible physical ori gin of this zone of avoid- 



Verde, Oh fc Jimenez (2002)) 



at tittp : //www. physics . rutgers . edu/ ^raulj/rot/roti . ps| 

(rotipseu.ps) 



The top-right panel of fig. [j] shows that low mass halos 
have on average higher concentration parameters than mas- 
sive halos, in agreement with theoretical predictions. This is 
not surprising since in hierarchical CDM models low mass 
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halos form at higher redshift than massive halos. The solid 
and dashed lines in the top right panel of fig. [I] show the 
median M — c relation predicted from numerical simulations 



(Bullock et al. (2001) figs. 4 and 5) for low density and 
high de nsity environments respectively. Although the recov- 



portion of the plot may be unstable to bulge formation and 
become bulge-dominated; disks may form preferentially out 
of low angular momentum gas, which settles in the center 
first; finally for some combinations of the parameters disks 
might fail to trigger star formation. In Verde, Oh & Jimenez 



ered distribution of c vs. M200 seems to show a somewhat (2002| ) we discuss in more detail the origin of the zone of 



steeper decline of c with increasing M200 (in agreement wit h 
the findings of |v an den Bosch, Burkert fc Swaters (2001 ) 



van den Bosch & l3waters™^2001 ) for 14 dwarf galaxies), it 



is difficult to draw a definitive conclusion due to the scat- 
ter in the theory prediction (simulations show a 30% scatter 
around this median relation) and in the recovered correla- 
tion. The overall distribution of halo concentrations is in 
reasonable agreement with ACDM predictions. This is more 
clearly seen in the left pa nel of fig. where the theoretical 
predictions by Jing (200C) are overplotted as a thick line on 
the distribution of c (obtained from those galaxies for which 
a meaningful fit was found). Maybe more interesting is the 
fact that very few galaxies are found to have values of c be- 
low the theoretical predictions. The small tail of 20 galaxies 
with c < 0.3 is mostly due to 14 galaxies for which fd > 0.2, 
i.e., above the value allowed by nucleosynthesis. For these 
galaxies the best fit parameter combination is thus unphysi- 
cal. In §4 we will argue that this might be due to limitations 
of our model for the rotation curves. 

The bottom-left panel shows that there is no evidence 
for a correlation between fd and c. On the other hand 
(bottom-right panel) there is a clear correlation be tween A' 



and fd- This is in agreement with that found by Burkert 
(2000|Qthc solid line shows his best least square fi t to the 
data) and |van den Bosch, Burkert fc Swaters (2001 ) but for 
only 14 dwarf galaxies. This is puzzling because cosmologi- 
cal models do not predict any correlation between A (where 
X = id/fdX) and fa . 

Burkert (2000) points out that this correlation can arise 
because of a correlation between the specific disk angular 
momentum (jd/ fd) and fd- If baryons that make up the 
disk cooled from the inside out, the outermost rings of the 
disk, that contain most angular momentum, cooled later; 
thus predicting a correlation with a slope very close to that 
observed. In this scenario the total specific angular momen- 
tum of the disk should be significantly smaller than that of 
the dark matter, since most of the galaxies have fd below 
nucleosynthesis value, but this seem s not to be the case for 
the 14 galaxies used in the ana lysis ( van den Bosch, Burkert 



& Swaters 2001 



Burkert 200 



On the right panel of fig. H we show the distribution of 
A' for galaxies of our larger sample for which a meaningful 
fit was found. Als o plotted as a solid thick line is the theo- 
retical prediction (Bullock et al. 2001) for the distribution of 
A found in N-body simulatio ns (which is very similar t o that 
predicted by linear theory ( Heavens fc Peacock 198^ )). Al- 



though the agreement is reasonable, there is a tail of galaxies 
for which the recov ered A' is below the predictions. 

Burkert (2000) did not find this, but t his may be due 



to the fact that the galaxi es in the sample of van den Bosch 



Burkert fc Swaters (2001) are LSB' s so they are selected to 
have high spin (Jimenez et al. 1998). 

We argue that the most significant feature of the 
bottom-right panel of figure 1 is the lack of galaxies above 
the solid line (what we call zone of avoidance). Many differ- 
ent effects might be at play here: disks in the bottom right 



avoidance. 

In Appendix II, we explore the degeneracies between 
model parameters. Although we find that there should be 
no significant degeneracies between model parameters for 
most of these high-quality rotation curves, we find that for 
the lowest quality rotation curves, significant degeneracies 
develop. The degeneracies have almost exactly the same di- 
rection as the 3 correlations recovered from the data. We 
therefore caution the reader that the correlations between 
model parameters seen above maybe the result of parame- 
ter degeneracies if there are systematic errors in the rotation 
curves or the quoted errors have been underestimated. 



3.2 Fits to disk models within a 
Pseudoisothermal halo 

We fit the pseudoisothermal halo model of §2.1 to the same 
set of galaxies as above using the same inverse variance 
weighted \ 2 . Our aim is two-fold: first to see if the same 
correlations among the parameters found for the NFW pro- 
file still hold for this new profile and second, to investigate 
the distribution of core sizes. 

We find that, for this dark matter profile model, for 114 
galaxies out of the 362 (i.e. 32%) the % 2 minimization proce- 
dure does not find a global minimum within the boundaries. 

Figure ^ shows the location of galaxies for the projec- 
tions of the three dimensional parameter space (M200, fd, 
c), where c = R^oo/fc- As for the NFW profile, there is a 
clear anti-correlation between fd and Mgoo and t he same 
zone of avoidance. In Verde, Oh & Jimenez (2002) we dis- 



cuss in detail the physical origin of this zone of avoidance. 
Here we only point out that this zone of avoidance seems to 
be independent of the shape of the dark matter profile. 

There is no evidence for a correlation between fd and c. 
As for the NFW profile we can now use the best fitting val- 
ues of (M200, fd and c) to compute the value of A' using eqs. 
19 and 20 from Appendix I. There is a clear correlation be- 
tween fd and A' (bottom right panel of fig. ^|), very similar to 
that found for the NFW model (fig. [!]). The distribution of 
A' for the whole sample is shown on the right panel of fig. ^| 
the thick solid line is the theoretically predicted value. Al- 
though the theoretical prediction is obtained from numerical 
simulations, where the dark matter halo profile is closer to 
the NFW, there is reasonable agreement. This is not sur- 
prising since the N-body predictions are very similar to the 
linear theory ones, and if halos were isolated (i.e., no further 
merging after formation) th ey would be well ap proximated 
by isothermal spheres (e.g.. 



Gunn & Gott (1972)) 



The upper-right panel shows that most galaxies in the 
sample have c = -R200 /r c larger than 20, already pointing to 
small sizes of the cores. This is better illustrated in the left 
panel of fig. where we show the histogram of R2oo/r c for 
the galaxies for which meaningful fit are obtained. There is 
a sharp cut-off at ifeoo/fc ~ 16, with 70% of the galaxies 
having i?20d/Vc > 100. This implies that the size of dark 
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0.0 0.5 1.0 1.5 2.0 2.5 3 
logio(R 2 oo/r c ) 



Figure 4. Pseudo-isothermal model. Left panel: distribution 
of the concentration parameter c = R20o/r c for the pseudo- 
isothermal dark halo. Note the sharp cut-off at c = 16, 70% 
galaxies in the sample have c > 100, i.e. 70% of galaxies have 
cores with sizes below 1% of the dark halo virial radius. Right 
panel: the distribution of A'. The continuous line is the numerical 
simulations prediction (see text for more details). 
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Figure 5. The size of the core radius in the pseudo-isothermal 
model compared to the break radius in the NFW model. Note 
that the former is 10 times smaller than the latter and that 70% 
of the sample has cores below 2 kpc. 



matter cores is modest, about 6% of the virial radius or 
smaller. 



4 GENERAL CONSIDERATIONS 

From the recovered best fitting disk parameters in the con- 
text of the 2 models (NFW profile and pseudo-isothermal 
profile) we can draw some general conclusions. First of all, 
the correlations between c and M200, fd and M200 and fd 
and A' are quite similar for the two dark matter profiles: 



these correlations and the zone of avoidance in the different 
projections of the three dimensional parameter space (M200, 
c, fd) seem to be independent from the shape of the dark 
matter profile. The disk parameter correlations we find from 
this large (362 galaxies) and diverse sample of galaxies are 
in agreement with those found from previous studies that in- 
volved m uch smaller samples (~ 20) and only dwarf or LSBs 



galax ies (Burkert 2000; van den Bosch, Burkert & S waters 



2001] ). 

Both NFW and pseudoisothermal profiles are good fits 
to about 70% the rotation curves. For the NFW profile we 
find no meaningful fit for 34% of the galaxies (122 out of 
362), while for the pseudoisothermal profile we find no mean- 
ingful fit for 32% (114 out of 362). We do not expect to find 
that the model provides a good fit for all the galaxies in 
the sample: galaxies might be more complicated than our 
model for the rotation curves. In fact some rotation curves 
show strong evidence for spiral arms and bars in the rotation 
pattern, and we have not attempted to model any bulge or 
bar component. However, we find that many rotation curves 
that have no meaningful fit for one dark matter profile are 
well fitted by the other: for only 52 galaxies (14% of the 
sample) there is no meaningful fit in either model. These 
correspond to the poorest quality rotation curves (low spa- 
tial resolution and poor spatial coverage, i.e. no flat or inner 
rotation curve). 

Fig. [5] shows that the core radius in the pseudo- 
isothermal profile (r c rso) is much smaller than the break 
radius in the NFW profile (r C NPw): fciso is about 10 times 
smaller than r C NPW- In other words, best fitting models for 
the dark matter in this sample of galaxies are those with 
steep dark matter profiles (r _1 or steeper), at least down to 
6% of the virial radius. 

The above considerations lead us to conclude that large 
dark matter cores (and therefore low central dark matter 
densities) seem to be excluded. In fig. ^ we show the distri- 
bution of central dark matter densities recovered from the 
best fits to the rotation curves. For the NFW profile model, 
since the central density is not finite, the densities have been 
evaluated at r = 0.5 kpc. This is comparable to the smallest 
radius at which the rotation curve can be reliably measured 
in the samples used in this papers. 

The ar row shows the central de nsity (at r = 0.5 kpc) 
inferred by Borriello fc Salucci (2001 ) (see their fig. 5) from 
9 galaxies. 

We hav e also selected 37 gal axies (out of the 48 LSBs) 
analyzed by ie Blok et al. (2001), that show indication of a 



flat core, and computed the distribution of inferred central 
densities (the 9 galaxies which show no flattening of the pro- 
file should have higher central densities). This distribution 
is shown as the dotted line in fig. ^ 

There are a few features worth commenting on. First, 
the mean density for the pseudo- isothermal profile is slightly 
lower than that for the NFW profile at 0.5 kpc. The aver- 
age value of deriv ed central densities fo r the 27 core galax- 
ies selected from ie Blok et al. (2001) is not significantly 



lower than the average value for the pseudo-isothermal pro- 
file. Therefore, there is quite a reasonable agreement for the 
central densities derived from our model and those derived 
inverting the Poisson equation in the above LSB galaxies. Of 
course, the NFW densities will be much higher at radii be- 
low 0.5 kpc. Since densities for smaller radii have not been 
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Figure 6. The distribution of central dark matter density for the 
two profiles at a distance from the center of 0.5 kpc. The dashed 
line corresponds to pseudo-isothermal dark halos while the solid 
line corresponds to NFW dark halos. Note the modest density of 
dark matte r. The arrow shows the c entral density (at r = 0.5) 
inferred by Borricllo fc Salucci (2001 ) and the dotted line is the 



distribu 
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ion of central densities for 27 galaxies from ie Blok et al 



t lat show a flat core. 



Figure 7. NFW model. The predicted (M/L) in our analysis, 
using the NFW dark halo profile, for the Courteau and PW sam- 
ples. The left panel corresponds to the Courteau sample. The right 
panel shows the histogram for the galaxies in the PW sample as 
a solid line. The dashed line shows the resulting distr ibution for 



the (M/L) rat io assuming a maximum disk model (Palunas & 



Williams 200C). The gray areas correspond to the allowed values 



from synthetic stellar populations. 



directly derived in a reliable way, it is difficult to assess 
whether or not the high dark matter densities predicted by 
the NFW profile are excluded. We simply point out that 
physical mechanisms that remove efficiently most of the dark 
matter inside the central kpc from a steep dark matter pro- 
file are currently under investigation (Merritt et al. 2002, in 
preparation; Milosavljevic et al. 2002, in preparation). 



4.1 M/L of the sample 

A good independent test of the values of the best fit disk 
parameters (in particular of fd) recovered from modelling of 
the rotation curves is a comparison of the baryonic mass- 
to- light ratios (M/L) inferred from the recovered disk mass 
M — fdM2oo and measured luminosity (L; from photom- 
etry) with values predicted by synthetic stellar population 
models. Stellar synthesis models assume an initial mass func- 
tion (IMF) and a star formation rate, age and metallicity 
to uniquely predict the M/L. The IMF is an ansatz which 
can be guided by observations in the solar neighborhood 
(e.g., Kroupa (2001)), while the star formation rate, age and 
metallicity can be inferred with help of the spectrum of the 
stellar population. In our case, because spectra are not avail- 
able, we have to make reasonable guesses for these param- 
eters. Once they have been fixed, current stellar synthesis 
models agree, quite remarkably, in their pr edictions about 
the M/L ratio (e.g., iBell & de Jong (2001b). However, un- 



certainty in the assumptions ultimately limit the precision 
in the M/L predictions from synthetic stellar population 



models leaving a 20-40 % uncertainty. Taking these uncer- 



t aint ies into account, a large disagreement between the M/L 
recovered from the rotation curve modelling and permitted 
values from stellar synthesis models, would imply that the 
modelling of the rotation curves is probably incorrect. 

Fig. M shows the baryonic M/L ratios derived by using 
the best fitting model values for M200 and fd for the NFW 
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Figure 8. Pseudo-isothermal model. The predicted (M/L) in our 
analysis and in synthetic stellar populations. The line styles and 
shaded areas have the same meaning as in fig. [?j 



profile (for the 240 galaxies for which meaningful fits were 
found) and the published values for the reddening-corrected 
i ntegrated photometry for the Courteau and PW samples 



ourteau 1997 



Palunas fc Williams 2000| ) ; Fig g shows the 



distribution of M/L ratios recovered by using the pseudo- 
isothermal profile model. In both cases, the recovered distri- 
bution of M/L is fairly broad. 

In both figures, the allowed range of M/L u sing the 
synthetic stellar population models developed in Jimenez 



et al. (1998) are shown as gray areas. This range is ob- 
tained by varying B — R (between and 1.5), the IMF 
(Salpeter; Miller-Scalo) and the star formation history of 
galaxies (a star formation recipe with a star formation rate 
M» tx exp(-t/t„), where 3 Gyr < U < 20Gyr). These as- 
sumptions allow 0.1 < (M/L) RI < 3 where the subscript 
R, I refers to the M/L in the R and I bands; the assumed 
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age of the stellar population has the greatest impact on the 
derived M/L range, with low values of M/L associated with 
progres sively younger stellar populations. 



We do not have extended information about the environ- 
ment of the galaxies from our sample, but very low valu es of 
fd are found also from analysis of LS B's galaxies (e.g 



Th 3 M/L derived from the NFW profile show a some- 



den Bosch, Burkert & Swaters (2001)) that are not found 



what broader range than allowed: for the Courteau sample 
4% of the galaxies have uncomfortably low M/L while 16% 
have M/L larger than allowed from synthetic stellar pop- 
ulations; for the PW sample 14% have low M/L and 7% 
have too high M/L. For the pseudoisothermal profile the 
situation is very similar: for the the Courteau sample 2% of 
galaxies have low M/L and 17% have too high M/L while 
for the PW sample 5% have too low M/ L and 10% too high. 

It is commonly claimed in the literature that NFW fits 
require M/L ratios which are too low; this is then used as 
an argument against dark matter p rofiles with central cusps 
(Ee Blok, McGaugh fe Rubin 20oij). The majority of studies 



have however focused on LSBs; we find that for HSBs only 
a small fraction of galaxies have excessively low M/L. Low 
values of M/L may plausibly be associated with very young 
stellar populations (e.g., if the last burst of star formation 
is ~ 1 Gyr old, then (M/L)n,i ~ 0.1. Tighter constraints 
on M/L are only possible with multi-color photometry or 
ideally, spectroscopy). Note that the distributions of M/L 
ratios derived from the pseudoisothermal profile are very 
similar to those derived using the NFW profile (fig. [j] and 8f) ; 
removing the central dark matter cusp does not change M/L 
significantly. Thus we conclude that the derived disk mass is 
fairly independent of the dark matter profile we assume and 
in broad agreement with values allowed by synthetic stellar 
population models. 

A possible criticism of studying HSB galaxies is that, 
in contrast to LSBs, the rotation curves of these galaxies 
are not dark matter dominated. In an extreme case, this 
would correspond to a maximum disk model, in which the 
circular velocity is attributable solely to the disk. In this 
case, we would be completely insensitive to the dark matter 
profile. The dashed line histograms in the right p anels of 



figures [j| and [] show the M/L ratio derived by 



Paluna; 



& Willams (200C) fitting a maximum disk model to their 
sample. Not surprisingly, their derived average M/L is larger 
than that of both the NFW and pseudoisothermal profile; 
however, it is also larger than allowed by the synthet ic stellar 
popula tion models . The Sa-Sd disk galaxies in the Palunas 
& Will jams (2000[ ) sample have B - R < 1.2, which fairly 
robustly implies (M/L)r i i < 2.5. The M/L ratios implied 
by the maximum disk models are much larger: they are more 
typically associated with ellipticals. 

4.2 Low recovered disk mass fractions fd 

It is worth commenting on the low values of fd we ob- 
tain. For both the NFW and pseudo-isothermal fits, we 



often obtain values of fd as low as fd 



even for 



M200 ~ few x 10 A/q. This is implausibly low: it implies 
only ~ 1% of the baryons in the halo cool to form the disk. It 
is difficult to think of any feedback mechanism (e.g., super- 
novae explosions) which could have such an extreme effect 
in inhibiting cool ing or expelling baryons from the disk (e.g., 
Efstathiou (200C )). A possible mechanism that can consid- 
erably reduce the amount of gas in disk galaxie s is RAM 
pressure stripping ( Quilis, Moore fc Bower 2000] ). However 



orbiting clusters. Thus the ram pressure stripping mecha- 
nism seems not to be responsible for the low values of fd we 
recover. 

These very low values of fd are also responsible for the 
tail of low M/L galaxies in figures [?] and ^. 

At present, we regard these low values of fd as pointing 
toward some inadequacy of the simple disk-halo model, or 
the assumed form of the dark matter profile. On the other 
hand, the high mass end of the galaxy mass function is dom- 
inated by elliptical galaxies and it is not inconceivable that 
we are selecting only those dark halos that manage to form 
a disk. 

Since the value of fd mainly affects the inner part of the 
rotation curve, we may ask: how should the dark matter pro- 
file be modified so that the recovered fd is more meaningful? 
It is clear that dark matter density in the center needs to 
be reduced. However, as discussed above, the best fit model 
for the pseudoisothermal profile has very small cores (and a 
steeper -r~ 2 - slope outside). 

We are investigating possible physical mechanisms 
within the CDM paradigm that can achieve such result, thus 
decreasing substantially the dark matter content of a NFW 
halo within the inner kpc (Merritt et al. 2002, in prepara- 
tion; Milosavljevic et al. 2002, in preparation). 



5 CONCLUSIONS 

We have analyzed a large set of high spatial resolution ro- 
tation curves of galaxies with the goal of determining the 
model parameters for a disk embedded within a cold dark 
matter (CDM) halo that we have modeled in two ways: ei- 
ther with a NFW profile or a pseudo-isothermal profile. To 
this aim we have developed the expression for the rotation 
curve of an exponential disk embedded in a pseudoisother- 
mal halo (Appendix I). 

Our study differs from previous work in two respects. 
Firstly, we use optical rotation curves rather than HI ro- 
tation curves, so the rotation curves we use are free from 
beam smearing and have smaller errors and higher spatial 
resolution. Secondly, while many studies have focused on rel- 
atively small samples of dwarfs, we use a large sample (400 
galaxies) spanning a wide range in luminosity and surface 
brightness. 

We find that the NFW profile provides a good fit to 66% 
of the g alaxies in the sample (in agreem ent with previous 
studies, van den Bosch & Swaters (2001)) that were based 



this mechanism works only on galaxies orbiting in clusters. 



on few (20) dwarf galaxies), with a distribution of recovered 
concentration parameters and spin parameters broadly con- 
sistent with that predicted by CDM numerical simulations. 
When the sample is fitted with an isothermal profile with a 
core, 68% of the galaxies are well fitted within this model 
and the best fit model favors cores with small sizes (almost 
all have r c < O.O6-R200 and ~ 70% have r c < O.Olifeoo)- 
However, we find that many rotation curves that have no 
meaningful fit for one dark matter profile are well fitted by 
the other: for only 52 galaxies (14% of the sample) there is 
no meaningful fit in either model. These largely comprise low 
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quality rotation curves. Our findings are consistent with the 
inner dark matter profiles being steep (slope of -1 or steeper) 
down to radii that are few percent of the virial radius. Large 
dark matter cores (and therefore low central dark matter 
densities) seem thus to be excluded. 

Numerical SPH simulations of disk formation within 
dark halos find that disks loose a significant fraction of their 
angular momentum if: 1) the dark matter profile is steep 
in the c enter and 2) the baryons sett le in right after virial- 
ization (Navarro & Steinmetz 2000a). On the other hand, 
we recover steep dark matter profiles and no significant loss 
of angular momentum. There are several routes to explain 
this. Some authors have pointed out that if baryons ar e al- 
lowed t o set tle well after virialization of the halo (e.g. 3aiz 
et al. (2001)), then there is no significant loss of angular 



momentum. Another possibility to consider is that current 
hydrodynamical simulations lack the resolution to simulate 
a multi-phase interstellar medium and also to follow the de- 
tailed formation of giant molecular clouds, and therefore the 
sites of star formation, within the settling disk. It is therefore 
not clear that current simulations contain all the physical in- 
gredients needed to simulate the complex process of baryons 
cooling down into a disk. 

In both models, for the galaxies that have meaningful 
fit to the rotation curves, the recovered baryonic mass to 
light ratios are broadly in agreement with predictions from 
synthetic stellar populations. 

From this large sample of galaxies we find several corre- 
lations among the best fitting model parameters and some 
regions in the parameters space that are not populated, we 
refer to those regions as zone of avoidance. These correla- 
tions an zone of avoidance are remarkably similar in both 
profiles and do not change when the rotation curves are fit- 
ted by a model with three free parameters or with four free 
parameters. In particular we find a strong correlation be- 
tween the disk mass fraction and the spin parameter for 
different kind of galaxies (HSBs and LSBs), in good agree- 
ment with findings of previous studies that involved only few 
dwarf LSBs galaxies. These correlations may be a result of 
parameter degeneracies if errors in the rotation curve have 
been underestimated; otherwise, they could have a physical 
origin. 

When fitting the sample with the NFW profile we ob- 
tain that for the 66% of galaxies for which meaningful fits 
to the rotation curves are found, the distribution of concen- 
trartion parameters is in good agreement with predictions 
from N-body simulations. Also, the distribution of disk spin 
parameters is in broad agreement with the distribution of 
spin parameters for the dark matter predicted from N-body 
simulations and linear theory, but with a small tail at low 
spins. 

The primary motivation in previous work of studying 
the rotation curves of dwarfs and/or LSB's was that they 
would be dark matter dominated. Although it is quite plau- 
sible, this claim must be taken with caution. In fact, our 
analysis of HSB and more luminous galaxies, as we have 
seen, yields values for the baryonic disk mass fraction very 
similar to previous analysis of dwarfs; a naive interpreta- 
tion of rotation curve fits alone would lead us to claim these 
galaxies as dark-matter domin ated as well. In Appendix II 



cannot be fit by NFW profiles, ~ 57% as opposed to ~ 33% 
for the entire sample (those which do fit have somewhat 
low but theoretically plausible concentration parameters c). 
This difference may be a 1 a statistical fluctuation (there 
are only 26 LSBs in the sample), it may reflect the greater 
fidelity of LSB rotation curves to the true underlying dark 
matter distribution, or it may perhaps be due to a system- 
atic difference in the dark matter profile of halos which un- 
dergo greater tidal torquing. On the other hand, for the 
pseudoisothermal model, only one LSB has fd > 0.2. This 
seems to indicate that LSBs are indeed fitted better with 
a dark matter core. For the LSB sample, 9 galaxies have 
dark matter cores with sizes below 1.5 kpc and 16 below 
3 kpc. This is in fair agreement with our previous findings 
that the size of the dark matter core is modest. It is also 
in agreement with the dark matter core size inferred from 
inverti ng the Poisson equat ion using the observed rotation 
curve (de Blok et al. 2001). Ultimately, our sample choice 
involves a trade-off: we gain increased precision in the ro- 
tation curve (which breaks parameter degeneracies) and a 
much larger sample encompassing a much broader class of 
galaxies, but also acquire a somewhat larger uncertainty in 
the disk contribution to the rotation curve. 

In general, LSB galaxies tend to have higher values of 
A' for a given fd and lower values of c for a given mass than 
HSB galaxies (see Appendix II). This supports the view that 
the dark matter profile of LSBs might reflect an enviroment 
in which dark halos undergo a higher tidal torque. 

If the M/L of the disk was well constrained, this ambi- 
guity can be eliminated: it would be possible to subtract the 
disk contribution to the rotation curve and recover directly 
the inner profi le of the dark halo b y inverting the Poisson 
equation (e.g., ie Blok et al. (2001)). The M/L can be con- 



we analyze the LSB sample of McGaugh, Rubin & de Blok 



(200l[[]and find that a somewhat larger fraction of them 



strained from multi-color photometry or even better, high 
S/N spectra of the stellar population in galaxies (extending 
as close as possible to the K band) together with high-spatial 
resolution rotation curves especially of the inner 1-2 Kpc. A 
more detailed treatment of this issue will be presented in a 
forthcoming paper. 
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APPENDIX I: ROTATION CURVE FOR 
PSEUDO-ISOTHERMAL PROFILE 

In this Appendix, we derive three relations used in the text 
for the pseudo-isothermal halo profile: (i) the rotation curve 
of the pseudo-isothermal sphere, equation (Q), (h) the rela- 
tion between the finite central densit y a nd the concentration 
parameter (c = r2oo/r c ), equation ( |15| ) ; (hi) a fitting for- 
mula for the disk scale-length Rd as a function of the halo 
mass and core-radius {M2oo,r c ), and the disk mass fraction 
and spin (fd, A') assuming angular momentum is conserved, 
equation (pc|). 

Starting from 



p(r) 



Po 



1 + (r/r c ) 2 
we obtain 



M(r) 



thus 



r' 2 p(r')dr = poAirr 2 



r c arctan 



(k) 



-arctan 



(k) 



and 

M200 = Po^nrl [c - 
where c = -R200 /r c 
M200 



arctan(c)] 

Since by definition 



A/3nR 3 200 
we have that 
200 

P0 = -— Pcrit — 

3 c - 



200p c 



arctan(c) 



(10) 



(11) 



(12) 



(13) 



(14) 



(15) 



The disk scale length Rd is obtained as follows. From 
the definition of the spin parameter A, the total angular 
momentum of the dark matter halo is: 



J : 



A 



\E\ 1 / 2 G- 1 M, 



-5/2 • 



(16) 



where E is the total energy of the halo. The total angular 
momentum of the disk Jd is given by: 



J d = 2n / V(r)^(r)r 2 dr=jdJ 



(17) 



where E(r) = So exp[— r/Rd], and the central density is 
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So = f d M 20 o/(2nR 2 d ). The total energy E can be ob- 
tained through the virial theorem 2K + (7 = 0, where k 
denotes the kinetic energy K and U the potential energy: 
E = k + U = (7/2, thus, 



£ = 



1 



2 2 0/ -t 4 

-Pot 8Gr c 



f GM(r 



(/;■ 



47rr rfr 



: arctan 



(18) 



(3 



Using equations (|lrj)-(|l9|) we obtain the following equation: 



arctan(a^) 2 / 

y —!-x 2 exp(-a;F)da;(19) 



where A' = Xjd/fd and Y = r c /Rd- Equation ( |19| ) can easily 
be solved numerically for Y and thus i?<j. We give here a 
fitting formula that is good to <J 5%: 



(20) 



where q is defined in eq^l9[ a = 2.3, /3 = 0.4745, 7 = 0.85, 
S = 1.045. Equation (GO) gives an expression for the disk 
scale-length Rd for a disk embedded in a pseudo-isothermal 
sphere. 



APPENDIX II: 
MODELS 



FIT TO 4 PARAMETER 



In this Appendix we explore how a model with four free pa- 
rameters (Af2oo, fd, A and c), i.e. without fixing Rd from 
observations, fares against the three free parameters model 
used in the paper. This is interesting for two reasons: it al- 
lows us to estimate how much worse the parameter estima- 
tion becomes when adding an extra free parameter and al- 
lows us to assess the robustness of the zones of avoidance. As 
in the paper we study both a NFW and a pseudo-isothermal 
dark matter halo. We also include the LSB sample and show 
that the trends for parameter correlations for LSBs are the 
same as for the Courteau and PW samples (HSBs). 

In this case, for each galaxy rotation curve we explore 
the entire likelihood surface within the same limits as the 
three free parameter model and 0.001 < A' < 0.5. 

In fig. [9] and |l^ we show the position of each galaxy in all 
six projections of the four dimensional space (fa, M200, A', c), 
for the NFW and pseudo- isothermal model s respectively. Di - 
amonds correspond to galaxies from the Courteau (1997|) 
sample 



Palunas 



Williams 



triangles to galaxies from 
(2000j[W squares to the LSB galaxies from [McGaugh, Ru 
bin & c e Blok (2o6~ij). 



First we discuss the case of the NFW profile (fig. |9|). 
The same correlations and zone of avoidance are present as 
in the three free parameter model; for 144 galaxies out of the 
400 no fit was found within the boundaries. The upper-left 
panel of fig. ^ shows the same zone of avoidance (low values 
of f d and high values of A' are avoided) as before and also a 
similar correlation between fd and A' although not as clear 
as the one in fig. |l| 

The change in the resulting distribution of recovered 
parameters can be seen clearly in the plane c — M200 ■ Here a 
large group of galaxies have recovered values of c and M200 in 
excellent agreement with those found in fig. ^. On the other 



hand, a clump with unphysical values of c < 1 is clearly 
visible. They correspond to galaxies with unphysically high 
values of fd (see the corresponding clump in the M200 — 
fd plane). About 25% of the rotation curves favours this 
different solution. These are low spatial resolution rotation 
curves: their quality is comparable to that of f 583-4 shown 
in fig. 

Note that in all 3 cases, the degeneracies between pa- 
rameters run in almost exactly the same direction as the 
correlations recovered from the data. These correlations may 
therefore be the result of parameter degeneracies if errors in 
the rotation curve have been underestimated. 

For the four free parameter case considered here the 
recovered M/L are slightly broader than for the three free 
parameter case. This indicates that for this case we are still 
recovering correctly the mass of the disk. We have computed 
the M/L for the eight systems in the LSB sample with reli- 
able photometry. For these systems, M/L = 0.1 — 2.5 in the 
B band for both the NFW and pseudo-isothermal profile, 
in fair agreement with the predictions from synthetic stellar 
population models. 

The recovered parameters for the pseudo-isothermal 
dark matter halo are shown in fig. [l^. Also in this case 
the zones of avoidance are preserved and the correlations 
among the different parameters are remarkably similar to 
those obtained from the three free parameter model. In this 
case for 141 galaxies out of 400, no fit was found within the 
boundaries for the parameters. 

From fig. [jij it can be seen that LSB galaxies tend to 
have higher values of A' for a given fd and lower values of c 
for a given mass than HSB galaxies. A similar trend occurs 
for the NFW profile. 

Figure |ll] shows the recovered values of Rd in this case 
plotted against the measured values from the photometry 
for the two profiles. We see that the disk scale length is very 
poorly constrained from the rotation curves alone. In par- 
ticular, for the pseudo-isothermal case the estimated disk 
scale-length is systematically several times larger than the 
true disk scale length. This corresponds to inserting a con- 
stant surface density sheet over the extent of the rotation 
curve. 



APPENDIX III: DEGENERACIES IN THE 
MODEL 

Similar modelling to that performed here, but with HI ro- 
tation curves, has shown that there are serious model de- 
generacies (e.g.. 



van den Bosch, Burkert & Swaters (2001)) 



For example, it is co mmon to expect a degeneracy betwe en 
M200 and c (see e.g., ie Blok, McGaugh fe Rubin (2001 )). 
In this section we argue that the main reason for this is the 
poor spatial resolution of the HI curves and not the model 
itself. To demonstrate this we will explore the likelihood sur- 
faces of the 400 NFW fits (same conclusions derived in this 
section apply to the pseudoisothermal fits). We choose to 
do this using fits to 4 parameters to demonstrate that even 
with such a large number of free parameters degeneracies 
are removed for high-resolution rotation curves. Obviously, 
for the case of only 3 free parameters, possible degeneracies 
will be even more restricted. 

First, it is useful to look at how the model behaves under 
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Figure 9. NFW model, four parameters fit. In all six panels diamonds correspond to the Courteau, triangles to the PW and squares to 
the LSB sample. The correlations between A and fd (upper- left panel), M200 and (middle-left panel) are very similar to that obtained 
for the three parameter fit. Note the distribution of c vs. M200- 



changes of parameter values. To illustrate this we choose a 
fiducial model with (M 200 , A', f d , c)=(lx 10 11 M Q , 0.05, 0.01, 
10) and explore some of the possible degeneracies present in 
the rotation curve. To do this we keep two parameters fixed 
and allow the other two to change from their fiducial value. 



This results in six possible combinations. There will be de- 
generacy if in any of these combinations the two rotation 
curves are indistinguishable. Fig. |l^ shows the percentage 
difference of the rotation curve so obtained with the fiducial 
model. Note that if the rotation curve is not known to better 
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Figure 10. Same as fig. |s| but for the model with disks inside pseudo-isothermal dark halos. Note that the correlations between A and 
ffl and between M200 and are very similar to those found for the NFW model and to those found for the three parameters fit. Note 
in the A' — panel that LSBs have higher values of A' than HSBs for a given value of f^. Also, in the c — M200 panel, LSBs have lower 
values of c than HSBs for a given mass. 
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Figure 11. Comparison of the recovered scale-length of the 
disk, when using a four free parameter model, for the pseudo- 
isothermal model (top panel) and the NFW model (bottom panel) 
with the observed stellar scale-lengths. The disk scale-lengths are 
very poorly constrained from the rotation curve alone. 
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Figure 12. The percentage difference between two rotation 
curves for which two parameters are changed at the time while the 
other two are kept constant in a fiducial model with (M20O1 A', 
fd, c)=(l X 10 11 , 0.05, 0.01, 10). The different curves correspond 
to variations of 5, 10, 15, 20 and 25% in the model parameters. 
The two parameters that are changed are shown on the top of 
each panel. The fiducial model is shown as a solid horizontal line. 
If the rotation curve is known with less accuracy than the typical 
variation among parameters, these will be degenerated. 



than 5%, then there are clear degeneracies in the A' — fd, 
A' — A/200 and fd — M200 planes. On the other hand, in the 
other planes there are no degeneracies if the rotation curve is 
well-constrained below 2 disk scale- lengths. Therefore, if the 
rotation curve is known to better than 5%, with significant 
spatial coverage, degeneracies among the parameters should 
not arise. Of the various parameter combinations, the most 
degenerate pair is A' — f d - 

We now turn our attention to the actual data. We first 
choose a good example to illustrate how well de generacies 



can be broken. We choose galaxy 11846 from the Courteau 
(1997^~^ample. Fig. [l3| shows the 95.4% confidence contour 
for two parameters having been marginalized, i.e. integrated 
the likelihood (C = exp[— x 2 /2]), on the remaining parame- 
ters (note that these quoted confidence intervals should not 
be over-interpreted: the true (larger) uncertainty is likely to 
be dominated by systematic errors and inadequacies in the 
model). The confidence contours are indeed rather small, 
implying errors in the recovered parameters of only about 
20%. From a Fisher matrix analysis (see section 5.1) we find 
that about 20% of the rotation curves in our sample have 
the above quality, 60% have ~ 50% errors in the recovered 
parameters. Only 20% have much larger error bars, with 



much poorer constraints on the recovered parameters; one 
of these galaxies is f583-4. 

The case of f583-4 is extremely useful to illustrate the 
large degeneracies that arise due to the lack of spatial res- 
olution in the rotation curve. This has only 9 points with 
large error bars. Fig. |l4| shows the 68.5% confidence con- 
tours. These are large and broad, and the parameters are 
very poorly constrained. This should not come as a surprise 
since one can hardly expect to constrain well 4 parameters 
with only 9 points in view of the results from fig. [l2| We 
stress here that 80% of the galaxies in the sample have much 
better resolution rotation curves thus the recovered param- 
eters are much better constrained. 

We explore the dependence of parameters on the spatial 
resolution of the rotati on curve by investigating i n further 
detail e376g2 from the Palunas fe Williams (2000[ ) sample. 
We perform the following manipulations to e376g2: first we 
have removed all points beyond 10 kpc (no-flat), to produce 
a rotation curve without the flat part. Next, to the original 
rotation curve we remove the inner points below a radius of 
3 kpc, so there is no rising part of the rotation curve (no-r). 

We also select a random subsample, only 1/3 of the 
original points, from the rotation curve, thus degrading spa- 
tial resolution (low- res). Finally, we fit independently the 
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Figure 13. The 95.4% confidence contours for galaxy 11846 in 
the Courteau sample for a pair of parameters marginalizing over 
the other two. Note the small size of the confidence regions due 
to the good spatial sampling of the rotation curve. 
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Table 1. Values of recovered parameters from the rotation curve 
of e376g2 after having performed several manipulations to its ro- 
tation curve (see text). The mass in units of 10 12 Mq. 



approaching and receding parts of the rotation curve. This 
gives us a handle on the effect of asymmetries in the ro- 
tation curve. The values for the parameters are shown in 
table |l|. Also shown are the values for the full fit (all) with 
errors obtained from marginalizing over the other three pa- 
rameters. The variation in the recovered values of A' is very 
small, while variations of up to 100% occur for the other pa- 
rameters. The larger error in the recovered values happens 
for the low-res case, i.e. when the spatial resolution of the 
rotation curve is decreased. It is interesting to see how the 
different parameters are affected. The largest error for M200 
is for the low-res case, while the largest error for c occurs 
when only one part of the rotation curve is used. For fa the 
largest deviation occurs in the low-res case. 

The three rotation curves used in the above example 
are shown in fig. [15l 



Figure 14. Same as fig. [l3| but the galaxy f583-4 from the de 
Block ct al. (2001) sample which has a poor spatially resolved ro- 
tation curve. In this case we show the 68.5% confidence contours. 
Due to the poor spatial resolution the confidence contours are 
much larger in this case. It is also useful to see that the regions 
of degeneracy are not inside the avoidance regions. 



5.1 Fisher matrix analysis 

When the galaxy sample is large, it might be inefficient to 
explore the whole likelihood surface as we have done in this 
case and maybe more important, it might not be possible in 
the case of very large datasets. Therefore it is useful to use 
a fast algorithm to estimate the error and, especially, the 
possible degeneracies. As an alternative to the exact study 
presented above, and as an illustration for a method that 
can be suitable for very large data sets, we apply the Fisher 
matrix analysis to our data set. 

We performed a Fisher matrix analysis of the likelihood 
function £, where \ 2 = — 21n(£). Recall that for M param- 
eters, the marginal error on each of the parameters can be 
obtained, under the assumption of Gaussian likelihood, from 
the Fisher information matrix 



d 2 ln£ 

d6id6i 



(21) 



where i and j run from 1 to M (in our case M = 4) and the 
angle brackets denotes the averaging over many realizations. 
The marginal error is then 



(22) 



Since our best fit parameters are those that maximize the 
likelihood, a Taylor expansion of In C around Oo gives 
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Figure 15. Three examples of rotation curves. The one on the left panel is from the Courteau sample, the one in the middle is from the 
PW sample and the one on the right panel from the LSB sample. 



In£(0o + A0)~ln£o + - 



1 d 2 ln£ 



(23) 



The elements of the Fisher matrix can be estimated as: 



and 



— 2 [ln£(0 o + Af^e, 



(24) 

ln£(0 o - A^ei) - 21n£ ] . 

(25) 



J ~ 2A0iA0 j 

[ lnL{0 + A6 t ei + AOjej) + 
ln£(0 o - ABtei - AOjej) - 
ln£(6» - A0 i e l + A^e^) - 
In £(0 O + A6» l e l - AOjej) ] . 

Furthermore, since the Fisher matrix is symmetric it 
can be diagonalized (F = UAU T , where A is a diagonal 
matrix with elements the eigenvectors Ai, Am) and a new 
orthogonal set of variables built such that N = U t O, where 
O is the original set of variables (A, c, fd and M200 in our 
case). Thus each row of U T gives the amount by which the 
original variables are mixed and therefore how degeneracies 
arise among them. 

We have performed the above Fisher analysis on the 
400 rotation curves, both for the NFW and the pseudo- 
isothermal profile. Our conclusions regarding degeneracies 
are in agreement with the findings of the full analysis per- 
formed on galaxies 11846 and f583-4 and illustrated in fig- 
ures |l| and 0. 



